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Abstract. In the real world, many complex systems are represented not by single 
networks but rather by sets of interdependent ones. In these specific networks, 
nodes in one network mutually interact with nodes in other networks. This paper 
focuses on a simple representative case of two-layer networks (the so-called duplex 
networks) with unidirectional inter-layer couplings. That is, each node in one network 
depends on a counterpart in the other network. Accordingly, the former network 
is called the response layer and the latter network is the drive layer. Specifically, 
synchronization between each node in the drive layer and its counterpart in the response 
layer (counterpart synchronization, or CS) of this sort of duplex networks with delayed 
nodes and noise perturbation is investigated. Based on the LaSalle-type invariance 
principle, a control technique is proposed and a sufficient condition is developed for 
realizing counterpart synchronization of duplex networks.Furthermore, two corollaries 
are derived as special cases. In addition, node dynamics within each layer can be 
various and topologies of the two layers are not necessarily identical. Therefore, the 
proposed synchronization method can be applied to a wide range of multiplex networks. 
Numerical examples are provided to illustrate the feasibility and effectiveness of the 
results. 

Keywords: complex network, duplex, counterpart synchronization, stochastic 
perturbation 


1. Introduction 

Complex networks abound in almost every aspect of science and technology. Examples 
include the Internet, the World Wide Web, social networks, metabolic networks, 
food webs, networks of citations between papers, among many others [D 121 i. 
Synchronization is one of the most common phenomena in nature that interacting nodes 
can reach a coherent state, and it has been extensively investigated and discussed in the 
past two decades n El n m m m do]. For example, Pecora et al. [7] used the master 
stability function (MSF) approach to analyze the stability of the synchronous state in 
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coupled systems, Huang et al. [S] classified synchronization into five categories based 
on the MSF approach, Wang and Chen investigated synchronization in small-world 
networks [9] and scale-free HH] networks. 

Many literatures, including the above-mentioned ones, are primarily focused on 
synchronization within single networks that do not interact with other networks. 
However, many real-world networks often interact with and depend on each other. For 
example, people in a society interact with each other via their family relationship, 
friendship, or formal work-related acquaintanceship m- Countries in the global 
economic system also interact via various international relations. Transportation 
depends on air traffic networks, railway networks and road traffic networks. Obviously, 
in describing and dealing with such problems, the multiplex network representation 
would be more appropriate than the single network. Not surprisingly, multiplex networks 
have attracted enormous attention in the past few years in various helds of application. 
For example, Xiong et al. [I2] analysed the correlation between the information diffusion 
process and the opinion evolution process and found obvious interaction between the two 
processes. Liu et al. [13] investigated preferred degree networks and their interactions, 
and found dramatically different behaviors between two very similar networks. 

Counterpart synchronization describes the individuals in one network behave 
coherently with their counterparts in other associated networks, so it represents 
harmonious coexistence of nodes in multiplex networks. This sort of synchronization 
in a duplex network can also be taken as the so-called outer synchronization between 
two networks and has attracted wide attention. For example, Wu et al. investigated 
generalized outer synchronization between two different complex dynamical networks 
by employing nonlinear control [H] . In particular, this synchronization has been widely 
applied in topology identification of complex networks. To name just a few, Wu [15] 
and Zhao et al. [16] employed complete outer synchronization to identify topologies 
for weighted complex networks, Zhang et al. HU et al. adopted generalized outer 
synchronization to recover network structures. 

Meanwhile, time delays are unavoidable in complex networks due to finite 
information processing and propagation speeds. They extensively exist in the real world, 
such as communication networks, gene regulatory networks, and electrical power grids. 
Time delays greatly influence behaviors of dynamical systems. Many literatures are 
focused on synchronization and control of complex networks with coupling delay among 
different nodes [mun]. 

Noise is another important factor affecting behaviors of dynamical systems, as it 
is inevitable due to environmental disturbance and uncertainties. Generally, noise is 
harmful. However, the presence of noise sometimes plays a positive role [20], such as 
in inducing synchronization [21] and in facilitating topology identification of complex 
networks [221123]. 

Motivated by above discussions, we investigate CS of duplex networks with delayed 
nodes and noise perturbation. Based on the LaSalle-type invariance principle for 
stochastic differential delay equations, we design adaptive controllers to synchronize 
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nodes of the response layer to their counterparts in the drive layer, and put forward 
some sufficient conditions for guaranteeing CS. 

The rest of this paper is organized as follows. Modeling of duplex networks and 
some preliminaries are introduced in Section |2l Sufficient conditions for CS in duplex 
networks are presented in Section [3l In Section 01 two numerical examples are provided 
to illustrate the feasibility and effectiveness of our method. Finally, some conclusions 
are drawn in Section [5l 

Notation: Some necessary notations used throughout the paper are introduced, 
(or A"*^) denotes the transpose of a vector x (or a matrix A), || x II 2 is the Euclidean- 
norm of X, (g) represents the Kronecker product, R” is the n-dimensional real space. 
In G R"^” represents an identity matrix of order n, C"'[a, b] (a, 6 G R, a < b) represents 
the n—order continuously differentiable function space in [a,b]. 

2. Modeling and preliminaries 

Consider a duplex network consisting of N nodes in each layer, as shown in Fig. [H 
For convenience, take the upper layer as the drive layer and the lower layer which is 
dependent on signals from the drive layer as the response layer. We are concerned 



Figure 1. A duplex network with unidirectional inter-layer couplings. Each node 
in the upper layer is unidirectionally connected to a counterpart in the lower layer. 
The topologies of the two layers can be different, representing the individual sort of 
intra-layer interactions. 

about the impact of delayed nodes and noise caused by control input. Thus a drive 
layer consisting of N linearly coupled nodes is described by 

N 

d:s.i{t) = [ii{t, Xi(f), Xi(f - r(f))) -h ^ CijFxj(t)]df, i = 1, 2,..., N, (1) 

i=i 








Counterpart synchronization of duplex networks 4 

and the response layer with control input is given by 

N 

dyi{t) = [ii{t,yi{t),yi{t - T{t))) + '^dijTyj{t) + Ui{t)]dt 

i=i 

+ai{t,ei{t),ei{t-T{t)))dw{t), i = l,2,...,N. (2) 

Here, Xi(t) = (x^, ...,Xi„)’^ G and yi(t) = (ya, e are state vectors, Ui(t) 

is the control input for node i, fj : M+ x M” x M” —)■ M"" is a continuously differentiable 
function determining the dynamical behavior of node i, T = {aij)nxn G is the 

inner coupling matrix, and C = {cij)]sfxN ^ is the coupling conhguration matrix 

representing the coupling strength and the topological structure of network ([I]), with c^ 
being dehned as follows: if there is a link from node j to node i {i ^ j), Cij 7 ^ 0 ; otherwise, 
Cij = 0. The diagonal elements of matrix C is ca = — for i = 1 , 2 , 

D = {dij)NxN £ is the coupling conhguration matrix of network ([ 2 ]), which has 

the same meaning as that of C. r(f) denotes time delay of nodes, ej(f) = yi{t) — Xj(f). 
The noise term in network ([2]) is utilized to describe perturbation caused by the 
control input process inhuenced by environmental huctuations [23]. In particular, 
CTj : M"*" xR"' xM" —)■ is called the noise intensity matrix, w(t) = {wi{t), ...,Wmit))~^ 

is an m-dimensional Brownian motion dehned on a complete probability space (fl, iF, P) 
with a natural hltration {Pt}t>o- 

Throughout this paper, we make the following assumptions: 

(HI) The noise intensity function x,y) {i = l,2,...,iV) satishes the Lipschitz 
condition and there exists positive constants p, q such that 

trace{aja^) < px~’^x + gy"''y. (3) 

Moreover a(t, 0, 0) = 0. 

(H2) There exists a positive constant M such that 
||fi(t, Xi(t),Xi(t - T(t))) - fi(t, yi(t), yi(t - r(t)))|| 

< M[\\{xi{t) - yiit))f + ||(xi(t - r(t)) - yi{t - T{t)))f]C (4) 

(H3) r(f) is a differentiable function with 

0 < f{t) < p < 1. (5) 

Obviously, this assumption is ensured if the delay rif) is constant. 

Our purpose is to design proper controllers so that the noise-perturbed response 
layer ([ 2 ]) can reach CS with the drive layer ([T]). For this purpose, some necessary concepts 
and a lemma of stochastic diherential equations are presented. 

Consider the following n-dimensional stochastic diherential delay equation: 

dz{t) = (fit, z(t), z{t — T))dt + ipit, z{f), z{t — T))dw (6) 

on t > 0 with an initial value f E r, 0], R”), where r, 0], R”) represents 

the family of all J3)~nieasurable bounded C'([—r, 0], R"')—valued random variables, the 
measurable functions 0, (p : [0, -|-cxo] x R" x R” —)■ R"’ satisfy the locally Lipschitz 
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condition and the linear growth condition. It is known that Eq. (jH]) has a unique solution 
for any initial value f that is denoted by z(t, on t > —r. 

Let X denote the family of all non-negative functions V{t,z) on 

M_i_ X M"', which are continuously once differentiable in t and twice differentiable in z. 
For each V G x the diffusion operator CV associated to (E]) acting on 

X M",R+) is defined by 


where dV/d 


dV dV ^ 1 , i 

z = {dV/dzi, ...,dV/dzn),d‘^V/d'^z = {d‘^V/dzidzj 


}nxn- 


(7) 


Lemma 2.1 (A Lasalle-type invariance theorem for stochastic differential equations 
m)- Assume that both (fit, u, v) and (pit, u, v) are locally hounded in (u, v) while 
uniformly bounded in t. Assume also that there are functions V G C^’^(R+ x R"',R+), 
7 G L^(R+,R+), anduji,uj 2 G C(R",R+) such that 


CVit,u,\) < 'jit) — Ci;i(u) -|-a; 2 (v), V(f, u, v) G R+ x R” x R"', 
a;i(u) > a;2(u), Vu G R” 

and 

lim inf I/(t, u) = oo. 

||u||^-|-oo 0<t<oo 

Then Ker{uji—uj 2 ) 7^0 and for every initial value f G C^^([—r, 0], R"'), the solution 
z(t, 0 ^S\) has the following property: 

lim dist{z{t; f), Ker{cji — a; 2 )} = 0 a.s. 

£—>-cx) 

Moreover, if Ker{ui — (^> 2 ) = 0, then for every f G ([—r, 0], R"), lim z(t;^) = 0 a.s. 


3. Sufficient conditions for CS of duplex networks 

In this section, we will first give the definition on top of duplex networks. 

Definition 3.1 The duplex network formed by the drive layer m and the response layer 
^ is said to almost surely achieve CS if 

ei(t) = lim ||yj(t) - Xi(t)|| = 0, i = 1, 2,..., iV. (8) 

£—)-oo 

With the network models and the definition given previously, we arrive at the 
following main theorem. 

Theorem 3.1 Let (HI), (H2) and (H3) hold. The response layer can almost surely 
achieve CS with the drive layer m with the following control scheme: 

N 

u* W = hij{f)Tyj{f) - ^i{t)ei{f),i = 1, 2,..., N, 

i=i 


( 9 ) 
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gi{t) = ki\\ei{t)f,bij{t) = -ej {t)Tyj{t)C = ( 10 ) 

where /c* > 0 (i = 1,2, are arbitrary constants, bij(t),gi(t) {i,j = 1,2, are 

adaptive parameters updating with network dynamics. 

Proof. Since ej(t) = yi{t) — Xj(i:), dynamics of the synchronization error between 
counterparts in layers ([T]) and ([2)) can be written as follows: 

N 

deiit) = [fi(t, yiit), yiit - r(t))) - ^i{t, Xi(t), Xi(t - r(t))) + '^{dijTyj{t) - CijT:sij{t)) 

j=i 

+Ui{t)]dt + ai{t, ei{t), ei{t - r(t)))dw(t), i = 1, 2,..., iV. (11) 


Consider the following Lyapunov functional: 

N N N N 


i=l j=l 


ij Cij) 


Y. Cs.it)-gf 


2=1 


2=1 


M 


N 


Y^Bj(e)e,(e)dB, 


( 12 ) 


Jt-r(t) 1 ~ If 

where ^ is a sufficiently large positive constant to be determined. Thus the diffusion 
operator C dehned in ((Tj) onto the function V along with the error system ffTTl) is: 

N N 

CV = 2^ef(t)[h(yi(t),yi(t - r(t))) - h(xi(t),Xi(t - r(t))) + ^(di^Ty^lt) - CijT^jit)) 

i=l j=l 

N N N 

+ - 2 + dij - Cij)eJ{t)Tyj{t) 

j=l i=l j=l 


N 


M 


N 


+2 - g)ej{t)ei{t) + ^ ej{t)ei{t) 

^ ^ i=i 


2=1 


1 — fi 


N 


N 


YC(t - f(t))ei(t - T(t)) + '^trace(aJ<T,). 


(13) 


With the well-known inequality 2x^y < x^x -f y''"y and Assumption (H2), one 
obtains 


2e7 {t)%{t, yi{t), yi{t - r(t))) - fi(t, Xi(t), Xi(t - r(t)))] 

< e7 {t)ei{t) + M[{eJ{t)ei{t) + ei{t - r(t))’^ei(t - r(t)))]. 

Let e(t) = (e^(t), eJ(t ),..., ejf{t)y, then 

N N N N 

CV < (1 + M)^e7(t)ei(t) + M^e7(t-r(t))ei(t-r(t)) + 2^^Cije7rej(t) 

2=1 2=1 2=1 j = l 

-2 5^^e7(t)ei(2) + ^ e7(t)ei(f) _ ^ e7(t - r(t))e^(t - r(t)) 

i=l ^ i=l d' 

N 

+ trace{aj CTj) 

2=1 
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= (1 + M + ^ )e^{t)e{t) + Me^{t — T{t))e{t — r(t)) + 2e^{t)Pe{t) — 2ge^{t)e{t) 

i — /i 

- r(t))e(t - r(t)) + pe^{t)e{t) + qe^{t - r(t))e(t - r(t)) 

1 — yU 

M P'*' + P 

< (1 + M + ^ ^ ^ + 2Amax(- 2 -) ~ 2^)e'''(t)e(t) 

+ M(r(t) —- r(t)) + pe'^(t)e(t) + ge'^(t - r(t))e(t - r(t)), (14) 

1 — /X 

where P = C (8) T. 

From Assumption (H3), one has < 0, which results in 

CV <{1 + 2A^ax( ^ + P - 2^)e’^(t)e(t) + ge’^(t - r(t))e(t - r(t)) 

^ J. yx 

=-a;i(e(t))+a; 2 (e(t - r(t))). (15) 

Let 


)T 


g>g*^ {l + 2Xr, 


P, 2M - uM 

—^ + p + q), 




(16) 


one gets Ci;i(e) > Ci; 2 (e) for any e 7 ^ 0. Moreover, lim inf V = 00 . From Lemma 

||e||—>■+00 0 <t<oo 


12.11 one obtains lim e(f; .^) = 0 a.s. for any initial data f G ([—r, 0], M”). This means 

t—>-oo ' 

that CS of the duplex network ([T]) and ([ 2 ]) can be almost surely achieved for almost every 
initial data. This completes the proof. 


Remark 3.1 In the duplex, the drive layer m and the response layer ^ may have 
different topologies. In addition, the configuration matrices C and D are not necessarily 
symmetric or irreducible, which means that the intra-layer topologies can be undirected or 
directed, and they may also contain isolated nodes and disconnected clusters. Therefore, 
the control scheme can be applied to a wide range of duplex networks with unidirectional 
couplings. 


Remark 3.2 It is obvious that when CS between the two layers fjfi and is almost 
surely realized, one has ej(t) ^ 0 as t —?> 00 for i = 1,2, ...,N. Furthermore, it renders 
gfit) —)■ 0 and hifif) —?■ 0 for i,j = 1,2,...,N. This means that gfif) and bij(t) will 
almost surely become constant. 


Based on Theorem 13.11 one can easily derive the following corollaries. 

Corollary 3.1 Assume that (HI) (H2) and (H3) hold. If the two layers have identical 
configuration matrices (C = O), then the drive layer m and response layer ^ can 
almost surely reach counterpart synchronization through the following simplified adaptive 
control: ufit) = -gi{t)ei{f),gi{t) = ki\\ei{t)\\. 

Corollary 3.2 Assume that (HI) (H2) and (H3) hold. If there is no noise 
perturbation, the duplex network m and can reach counterpart synchronization 
through the following adaptive control: Uj(f) = Yl^=i^ij{t)Tyj{t) — gi(t)ei{t), gfit) = 
ki\\ei{t)\\‘^,bij{t) = -ejTyjit). 
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4. Numerical simulations 


In this section, two examples are given to illustrate the feasibility and effectiveness of 
the proposed synchronization scheme. 


Example 4.1 Consider a duplex network, each layer being composed of 5 nodes. The 
chaotic Lii system with various parameters is taken as node dynamics, with the i—th 
(i = 1,2, ...,5) node in both layers being described by 


±i = fi{t, Xi(f), Xi(f - r(f))) = 


f {36 + i*0.1){xi2{t) - Xii{t)) \ 

-Xii{t - T)xi3{t - r) + 20xi2{t) 



y Xifit - T)Xi2{t - t) 

- 3x^3 (f) / 

(36 + i * 0.1) 

36 + i * 0.1 0 ^ 


^ Xil ^ 


^ 0 

0 

20 0 


Xi2 

+ 

-Xii{t - T)Xi3{t - T 

0 

0 -3 y 


\ Xi3 y 


y Xii{t - T)Xi2{t - t) 


= Hxi{t) + G{xi{t - r)). 


(17) 

\ 

/ 

(18) 


Since the Lii system is chaotic, it is bounded in a certain region IWf . Thus there exists 
a positive constant R such that \\yk\\ < R and \\zk\\ < R for k = 1,2,3. Therefore, one 
has 


||G(y) - G(z)|| = \/[z3{yi - zi) + yiigj^ - Z3)f + [yi{y 2 - ^ 2 ) + ^ 2 ( 1/1 - ^ 1 )]^ 
< \/2i?||y - z||. 


(19) 


That is to say. Assumption (H2) is satisfied with M = y/2R for i = 1, 2,..., 5. 

The configuration matrices C and D for the drive layer and response layer are given 


as 


C = 


/ -6 
3 
0 
3 


V 1 


2 

-4 

1 

0 

0 


0 

1 

-4 

3 

0 


3 

0 

3 

-7 

2 


1 \ 
0 
0 
1 

-3/ 


and D 


/ -3 0 2 0 1 \ 

0-613 2 

3 1-40 0 

0 10-32 

\ 1 2 0 1-4/ 


( 20 ) 


respectively. The inner coupling matrix is taken as r=[l 1 0;0 1 0;0 0 1] and node 
delay is r = 0.003. Take afit^Ci^efit — r)) = aodiag^enit) — enit — r),ei 2 (t) — 
eait - T),ei 3 {t) - ei 3 {t - t)), (To = 1 for i = 1,2,..., 5, then ai{t,ei,ei{t - t)) 
satisfies the Lipschitz condition and the linear growth condition. That is, trace{ajat) < 
2aQef (t)ei{t) + 2aoef (t — T)ei(t — T). Meanwhile, assume that w(t) = [tci(t), tC 2 (t),'M^ 3 (t)] 
is a three-dimensional Brownian motion. The initial values of the i—th nodes in the drive 
and response layers are set to be {xnit), xait), Xioit)) = (1+ 0.3f,—0.6 + 0.3f, 0.3 + 0.3f), 
and {yii{t),yi 2 {t),yi 3 {t)) = (1 - sini, 1 - 0.3 cosi,-0.3i), f G [-r, 0] {i = 1,2,..., 5), 
respectively. The initial values of adaptive gains gfit) {i = 1,2, ...,5) and adaptive 
parameters bij(t){i,j = 1,2, ...,5) are chosen randomly in (0,1). 

Figure\^ shows the counterpart synchronization error of the duplex network (QP and 
The left panel shows Cijifi), while the right panel shows the total synchronization 
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error || e ||= y *■5 obvious that CS is almost surely 

achieved once the proposed control scheme is employed. Figure 0 further displays the 
adaptive feedback gains gi{t) {i = 1,2,..., 5) and adaptive parameters bij{t) {i,j = 
1,2,..., 5.) varying with time. It is seen that all the parameters reach constant values, 
which is consistent with Remark \3.S\ . 




Figure 2. The counterpart synchronization error between the drive layer m and the 
response layer ([2]) formed by Lii oscillators. Left: Cijit) varying with time t; right: the 
total synchronization error. 




Figure 3. The adaptive feedback gains gi{t) [i = 1,2,..., 5) (left) and parameters 
bij{t) {i,j = 1,2,..., 5) updating according to (ITUl) (right). 

Example 4.2 Synchronization in neuronal networks is one of the burning problems in 
neuroscience in recent years, and the Hindmarsh-Rose model l27^ has become a popular 
model for analysis of neuronal activity and has also been extensively investigated. For 
example, Fang et al. investigated chaotic synchronization of nearest-neighbor diffusive 
coupling Hindmarsh-Rose neural networks in noisy environments fMEf, and Zhou et 
al. discussed the Hindmarsh-Rose model by using impulsive pinning control fMSj. In 
what follows, we will discuss the synchronization between two coupled Hindmarsh-Rose 
neuronal networks. The Hindmarsh-Rose model can be described by a three-dimensional 
nonlinear differential eguations as follows HI- 

Xi = f(t, Xi(t), Xi(t - r(t))) 
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t 

Figure 4. The counterpart synchronization error of a duplex Hindmarsh-Rose neural 
network. Left: eij{t) varying with time t; right: the total synchronization error. 


/ Xi2{t) - Xisit) - Xii{t - t)^ + 3xii{t -rf + I ^ 

= 1 - Xi2{t) - 5xa{t - rf 

y ia{A{xii{t) +x) - Xi3{t)) ! 

Take / = 3,x = 1.56,// = 0.006, (Tq = l,r = 0.1. Assumption (H2) is satisfied 
The inner coupling matrix F = [1 1 0; 0 1 0; 0 0 1]. The intra-layer topologies, the 
noise term and initial states of nodes are taken as the same as those in the previous 
example. Figure\^shows counterpart synchronization errors between two unidirectionally 
connected Hindmarsh-Rose networks. Figure 0 further presents the updated feedback 
gains gfit) {i = 1,2,..., 5) and adaptive parameters bijifi) {i,j = 1,2, ...,5). It is clearly 
seen that the numerical simulations perfectly match the theoretical results. 




Figure 5. The adaptive feedback gains gfit) [i = 1,2,..., 5) (left) and parameters 
ihj = 1; 2,..., 5) updating according to (ITOl) (right). 


5. Conclusions 

In this paper, counterpart synchronization of duplex networks with delayed nodes 
and noise perturbation has been investigated. Based on the LaSalle-type invariance 
principle for stochastic differential equations, a sufficient condition guaranteeing CS 
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with the proposed control scheme has been provided. Numerical examples have also 
been presented to illustrate the effectiveness of method. The proposed method will find 
its applicability to a wide range of practical duplex networks. 
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